Cavity QED and Quantum Computation in the 
Weak Coupling Regime II : Complete Construction 
of the ControUed-ControUed NOT Gate 

in 

O : Kazuyuki FUJII t Kyoko HIGASHIDA } Ryosuke KATO } Yukako WADA 

^ ■ Department of Mathematical Sciences 

^ '. 

Yokohama City University 
Yokohama, 236-0027 

> ■ 

■ Japan 

O 

o 
o 

:^ 

^3 

> 

■ We solve the Schrodinger equation of the model in the weak coupling regime 
to construct the controlled NOT gate in the case of n=2, and to construct the 
controlled-controlled NOT gate in the case of n=3 by making use of several res- 
onance conditions and rotating wave approximation associated to them. We also 
present an idea to construct general quantum circuits. 



Abstract 

In this paper we treat a cavity QED quantum computation. Namely, we consider 
a model of quantum computation based on n atoms of laser-cooled and trapped 
linearly in a cavity and realize it as the n atoms Tavis-Cummings Hamiltonian 
interacting with n external (laser) fields. 
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The approach is more sophisticated than that of the paper [K. Fujii, Higashida, 
Kato and Wada, Cavity QED and Quantum Computation in the Weak Couphng 
Regime, J. Opt. B : Quantum Semiclass. Opt. 6 (2004), 502]. 

Our method is not heuristic but completely mathematical, and the significant 
feature is based on a consistent use of Rabi oscillations. 

1 Introduction 

Quantum Computation (or Computer) is a challenging task in this century for not only 
physicists but also mathematicians. Quantum Computation is in a usual understanding 
based on qubits which are based on two level systems (two energy levels or fundamental 
spins) of atoms, See as for general theory of two level systems. The essence of Quantum 
Computation is to construct an element of huge unitary group U{2"') by manipulating n 
atoms with photons, laser fields, etc. 

In a realistic image of Quantum Computer we need at least one hundred atoms. How- 
ever, we may meet a very severe problem called Decoherence which destroys a superpo- 
sition of quantum states in the process of unitary evolution of our system through some 
influence arising from natural environment. At the moment it is not easy to control the 
decoherence. See for example the papers in 2j as an introduction to the problem. 

An optical system like Cavity QED may have some advantage on this problem, so 
we want to consider a model of quantum computation based on n atoms of laser-cooled 
and trapped linearly in a cavity. As an approximate model we realize it as the n atoms 
Tavis-Cummings Hamiltonian interacting with n external (laser) flelds. As to the Tavis- 
Cummings model see fJ]. To perform the quantum computation we must flrst of all 
show that our system is universal 'S,. To show it we must construct the controlled NOT 
operator (gate) explicitly in the case of n = 2, [S], [H]. 

For that we must embed a system of two-qubits into a space of wave functions of the 
model and solve the Schrodinger equation. In a reduced system we can construct the 
controlled NOT by use of some resonance condition and the rotating wave approximation 
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associated to it. Then we need to assume that the couphng constants are small enough 
(the weak coupling regime in the title). 

Next we want to construct the controlled-controlled NOT operator in the case of 
n = 3. For that purpose the construction of controlled NOT gates of three types is 
required ^ because three atoms are trapped linearly in the cavity. We have given an idea 
to construct them in [3^, which is however not complete. In this paper we give a complete 
construction to the controlled-controlled NOT gate. 

However, to push on with our method (namely, for the case of n atoms (n > 5)) is not 
easy by some severe technical reason. Therefore we present the idea in once more to 
construct general quantum circuits, which will give a general quantum computation. 

The contents of this paper are as follows : 

Section 1 Introduction 

Section 2 A Model Based on Cavity QED 

Section 3 Quantum Computation 

3.1 Controlled NOT Gate 

3.2 Controlled-controlled NOT Gate 
Section 4 Further Problem 

Section 5 Discussion 



2 A Model Based on Cavity QED 

We consider a quantum computation model based on n atoms of laser-cooled and trapped 
linearly in a cavity and realize it as the n atoms Tavis-Cummings Hamiltonian interacting 
with n external (laser) fields. This is of course an approximate theory. In a more real- 
istic model we must add other dynamical variables such as positions of atoms and their 
momenta etc. However, since such a model is almost impossible to solve we consider a 
simple one. 

"'^In the study of Cavity QED Quantum Computation this important point is missed 
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Then the Hamiltonian is given by 

n 

5^ h, (aj+)e'(^^*+<^^) + aj-)e-'(^^*+<^^)) ® 1 (1) 

where is the frequency of radiation field, A the energy difference of two level atoms, a 
and are annihilation and creation operators of the field, and g a coupling constant, VLj 
the frequencies of external fields which are treated as classical fields, hj coupling constants, 
and L = 2". Here crj^'', crj and are given as 

aj'^ = I2 ® ■ ■ • ® I2 ® ® I2 ® ■ ■ • ® I2 (j - position) G M(L, C) (2) 
where s is +, — and 3 respectively and 

''-[\-\ ''-[V\ 

See the figure 1 as an image of the model. Note that the cases of = 2 and 3 are our 
target through this paper. Here we state our scenario of quantum computation. Each 
independent external field generates a unitary element of the corresponding qubit (atom) 
like ai (S> 02 ® ■ ■ ■ ® a„ where Oj G f/(2) , while a photon inserted generates an entanglement 
among such elements like c'-i ® 0-2 ® ■ • • a„. As a whole we obtain any element in f/(2") 
(a universality). 

Here let us rewrite the Hamiltonian (^. If we set 

n n 1 " 

s^ = Y.^';\ S- = Y.^^-\ s, = W^f\ (4) 

then can be written as 

H = ulL®a)a + AS3 (g) 1 + 5f (5*+ ® a + 5'_ (g) a"^) + 

n 

J2 hj (af ^e*(^^*+^^) + aj"^e-'(^^*+'^^)) 1 = Ho + \/(t), (5) 

which is relatively clear. Hq is the Tavis-Cummings Hamiltonian and we treat it as an 
unperturved one. We note that {5+, S^, S3} satisfy the su(2)-relation 

[^3, S+] = S+, [53, 5_] = 5_] = 2^3. (6) 
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Figure 1: The general setting for a quantum computation based on Cavity QED. The 
dotted hne means a single photon inserted in the cavity and all curves mean external 
(laser) fields (which are treated as classical ones) subjected to atoms 

However, the representation p defined by 

is a full representation of sm(2), which is of course not irreducible. 
We would like to solve the Schrodinger equation 

ZjU = HU={Ho + V)U, (7) 

where U is a unitary operator. As an equivalent form let us change to the interaction 
picture, which is performed by the method of constant variation. The equation 
i-^U = HqU is solved to be 



Un 



where = a^a is the number operator and Uq a constant unitary. Here we have used the 
resonance condition 

uj = A (8) 

, see for example jZj. By changing Uq i — > Uoit) and substituting into ((Z|) we obtain the 
equation 

after some algebras. This is the interaction picture of (|7j) and we use this for our quantum 
computation. Therefore we must calculate the right hand side of © explicitly, which is 



however a very hard task due to the (comphcated) term e ^^9{s+m+s (g)at) ^ convenience 
in the following we set 

A = An = S+ a + S_ (10) 

It has been done only for n = 1, 2, 3 and 4 as far as we know, see [7], [H]. We list the 
calculations for n = 1, 2, 3 in the following, 

One Atom Case In this case A in ()10|) is written as 

/ 



Ai 



By making use of the simple relation 




Bl/2 



'111 





' aa^ ^ 




( N + l 








- 


") 


1^ i^a j 




\ 


N ) 



(12) 



we have 



CO ( 1 \n oo ( -| \n 



{tgY-A^-rY. 



cos{tgVNTl) -.^^^ 



o(2^ + l)' 
a 



—i- 



.SmitgVTV) 



-a' cos 



(13) 



C{N + l) -iS{N + l)a 
-iS{N)a^ C{N) 

We obtained the explicit form of solution. However, this form is more or less well-known 
see for example the second book in [1 . 

Two Atoms Case In this case A in ()10|) is written as 

'^0 a a ^ 



Ao 



at a 
at a 
at at 



(14) 
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Our method is to reduce the 4 x 4-matrix A2 in ()14|) to a 3 x 3-matrix Bi in the 
following to make our calculation easier. For that aim we prepare the following matrix 






1 








1 

V2 





1 

V2 





1 

V2 





1 

V2 














1 



then it is easy to see 



/ 



T^AoT 



\ 



(15) 



V2a 
V2a^ V2a \ Bi 
\ V2at ^ 

where Bi = ® a + J- ® a) and { J+, J_} are just generators of (spin one) irreducible 
representation of Q- We note that this means a well-known decomposition of spin 



I ® I = © 1. From the decomposition we have 



-itgA2 _ '-p 



Tt. 



-itgBi 



(16) 



Therefore to calculate e ^^^^^ we have only to do e . Noting the relation 



and so 



Bl 



2(A^ + 1) 

2{2N+1) 
\^ 2(at)2 2N 



( 



2(2iV + 3) 



2(2Ar + 1) 



Bi = DBi 



2{2N -1) J 



(17) 



2n+l T-)n 
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we obtain by making use of the Taylor expansion 



CO ( 1 \n 



-itgBi = 1 + ^ 



1 \n oo / -| \n 



nt^,(2n+l)! 



-ih{N)a^ l + 2f{N) -ih{N)a 



where 



-1 + cos ft^./2(2A^ + 1)) sin (tgJ2(2N+l) 

f{N) = 1^, h{N) = - 



V2N+1 



Three Atoms Case In this case A in (jlOj) is written as 



A, 
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We would like to look for the explicit form of solution like (fTT^ . (fTK|l . If we 
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then it is not difficult to see 

^ a 

a) 



a 
at 



VSa 

VSat 2a 

2at VSa 

V3at 



1/2 



5 



3/2 / 



This means a decomposition of spin \ ®\ ®\ = 'k®i^®^- From the decomposition we 



have 



/ p-itgBi/2 



-itgAa _ rp 



\ 



V 



Tt. 



(21) 



^-itgBg^2 

Therefore we have only to calculate e"**^'^^/^ ^ which is however not easy. In this case there 
is no simple relation like ()12p or p7|) . so we must find another one. 
Let us state the key lemma for that. Noting 

^ 3N + 3 2V3a^ ^ 

7N + 4: 2v^a2 

2V3(at)2 7N + 3 

\^ 2\/3(at)2 3N 

VS{7N + ll)a 

20(iV + l)a 

20A^at 



-^3/2 



r3 

-°3/2 



y3(7A^ + 4)at 



6a^ 


V3{7N + 3)a 



V 



and the relations 







VS{7N - 4)a"f 







rj2n+l 75 D2n 752n+2 r}2 r>2n 

-°3/2 — -03/2-03/2; -D3/2 — -D3/2-D3/2' 
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we can obtain B^J^ and B'^J^'^ like 



B. 



2n 
3/2 



an{N + 2) 

Pn{N + l) 

2^UN-l){o)? 





ln{N) 





2v^en(iV+l)a2 



R2n+1 



v/3/5„(Ar + 2)a Qin{N + 2)a^^ 

^/?>(3n{N + l)a) 2en+i(A^+l)a 

2e„+i(iV)at y37n(iV)a 



1^ QUN-l){a^f 







t 







where 



a„(iV) = {v+\l - v.\l)/{2^), (3n{N) = {w+\l - u;_A^)/(2v^), 
7„(iV) = KA^-! - t;_A^'[)/(2v^), 5„(iV) = {w+\l - u-_A'^)/(2v^), 
en(iV) = (A!t-A!!)/(2v^), 

and A± = A±(A^), v± = v±{N), w± =w±{N), d = d{N) defined by 



5n(iV-l) / 



(22) 



(23) 



A±(iV) = 5A^ ± ^d{N), v±{N) = -2N - 3 ± ^d{N), w±{N) = 2N - 3 ± ^d{N), 
d{N) = IQN^ + 9. 



Then by making use of (j22j) and (jzHj) we have 



( 



OO f 1 ^''^ '-^ 



(2^)' 
/2(iV + 2) 



-■3/2 



o(2^ + l)' 



-V^iFi{N + 2)a 2V3hi{N + 2)a^ -6iHo{N + 2)c 
-2iHi{N + l)a 2y/3hi{N + l)a^ 



\ 



-y3zFi(iV + l)at /i(iV + l) 

2v^/ii(A^)(at)2 -2i/7i(A^)at fo{N) -VSiFo{N)a 

\ -QiHo{N - l){a^f 2VShi{N - l){a^y -V3iFo{N - l)a^ f_i{N-l) ) 

(24) 
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where 



/2(iV) 

fo{N) 
f-i{N) 
h{N) 

Fi{N) 

Fo(iV) 
HiiN) 
HoiN) 



v^{N)cos{tg^KiN)) - v^{N)cositg^X^{N))j /{2^d{N)), 
w4N)cositg^\4N)) - w4N)cositg^\4N))^ /i2^/d(N)), 
v^{N)cositg^\4N)) - v^iN)cositg^KiN))^ /i2^/d(N)), 
w+{N)cos{tg^\^{N)) - w4N)cos{tg^X+{N))] /{2^d(N)), 



cositg^X^N)) - cos{tg^X4N)) /{2^d{N)) 



^A_(iV) ^A+(iV) J 

^X+{N)sin{tg^X^{N)) - ^ X^{N)^m{tg^ X^{N))\ /(2v^), 



Now we must calculate the term 



F{t) = e^*3^" (e 



e**"'Mr(t) e 



from Q, so we introduce a brief notation 



V{t) ^ 



Therefore we want to calculate 



(25) 



-itgAn 



T(block diagonals)T^V^(t)T(block diagonals)^T^ 



explicitly for the case of n = 2 and 3. For that let us calculate V{t) and T'^V{t)T in 
advance. The calculation is straightforward and the result is 



Two Atoms Case 
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V{t) 



\ 



q(t) pit) 
q(t) p(t) 
Pit) q{t) 
\ Pit) qit) J 
with pit) = h^e^m+^)t+<t>i} ^ q^t) = /i2e*^(^^+'^)*+'^2> 



and 



T^Vit)T = 






-P+1 
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V2 





V2 






where we have omitted the time t for simphcity. 
Three Atoms Case 
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with pit) = h^e^^^^^+'^^t+M ^ qi^t) = /i2e^«"2+'^)*+'^2}^ ^ /ige^^"^ 
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and 



T^V{t)T 
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p+g+r 
V3 



2(p+g+f) 
3 



V6 ^ 
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2(p+g+r) 
3 



p+g+r 
V3 



n ^ 



2p—q—r 




p+g+r 
V3 



(29) 



where we have omitted the time t for simphcity. 



To write down all components of F{t) = e**^"^" \/(t)e~**^^" is very long and we moreover 
don't need all of them, so we omit it here. Next let us go to our quantum computation 
based on a few atoms of laser-cooled and trapped linearly in a cavity. 



3 Quantum Computation 

To develop a quantum computation based on atoms laser-cooled and trapped linearly 
in a cavity (Cavity QED Quantum Computation) a quick and clear construction of the 
both controlled NOT gate and controlled-controlled NOT gate is required, 5 . Let us 
construct such very important quantum gates in this section. 

3.1 Controlled NOT Gate 

In this subsection we treat the case of two atoms (the system of two qubits). First let us 
make a short review of the system of two-qubits. Each element can be written as 

ip = a++|+) 1+) + a+_|+) ® |-) + a_+|-) ® |+) + a_„|-) ® |-) 
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with two bases |+) and |— ) and |a++p + |a+_p + |a_+p + \a p = 1. Here if we identify 



1+) 



1 





then ip above becomes 



-) 



a I 



(30) 



How do we embed two-qubits in our quantized system ? It is not known at the current 
time, which may depend on methods of experimentahsts. Therefore let us consider the 
simplest one like 



® |0), 



(31) 



^ a++(t) ^ 
a+_(t) 
a_+(t) 

where |0) is the ground state of the radiation field (a|0) = 0). We note that in full theory 
we must consider the following superpositions 



\m) = E 

n=0 



® \n) 



«+-,n(i^) 

a-+,n(i) 

as a wave function, which is however too complicated to solve. 

To determine a dynamics that the coefficients will satisfy we sub- 

stitute (pTTjl into the equation 



zj^im) = m\m) = e^''^w{t)e-''^^^\m) 



1 



T^V{t)T 



1 



-itgBi 
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Let us rewrite the above equation. If we set 





( inJf] \ 






,r, ( + \ 


® |0) = 
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V2 



a++(t) 

a+_(t)+a_+(f) 

a__(t) J 



® |0) 



(32) 



then 



T^V{t)T ^ 



-itgBi 



(33) 



On the other hand, we have calculated the term T''V{t)T in ()27|) . 

Note that the above equation is not satisfied under the restrictive ansatz (jH^ . Because 
the left hand side contains only the ground state |0), while the right hand side contains 
the ground state |0) and some excited states |1), |2), |3). However, the states |1), |2), 
1 3) which have no corresponding kinetic terms contain the coupling constants hi and /12 
(see p{t) and q(t) in (j2EI)), so the equation is approximately satisfied if they are small 
enough (namely, in the weak coupling regime in the title). 

Therefore the (full) equation is reduced to the equations of {ipo, ipi, ip2, v^s} at the ground 
state |0) : 



.d 

I— 
dt 
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V Vsit) ) 



(34) 



where 



X12 — X21, 
p-q 



Xu 

X21 
X23 



V2 ' 
-p + q 

p + q 



1 + 3/(1: 



^ ^1 + 2/(0) + |/(1) + i/(0)/(l) + A(0)A(1) 



2^32 — 3^23; 
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and 



^34 = ^(1 + 2/(0)) 
X41 = Xi4, 
3^43 = 2^34 



-1 + cos{y/2gt) 



/(O) = 

h{0) = sin(v^^t), h{l) 



^^^^ ^ -1 + cos{y/Qgt) 
sm{^/6gt) 
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How do we solve it ? We use some resonance condition and the rotating wave approx- 
imation associated to it, which is popular in quantum optics or in a field of laser physics. 
Let us focus on the (2,3)-component of the matrix which came from the interaction term. 
The products /(0)/(l) and h(0)h(l) contain the term e~**3{v^+v^) Euler formulas 

cos(^) = (e^^ + e^'^)/2, sm{9) = (e*^ - e-'^)/2i. Noting 

^i{{ni+Lu)t+(f>i}^-Ug{V2+V6) ^ ^i{{ni+LU-{V2+V6)g)t+(t>i} 

we set a new resonance condition 

ni+uj-{V2 + V6)g = 0. (35) 

All terms in except for the constant one Q'{i^i+'^-i^+^)9)t+'t>i} = q^<I>i contain ones like 
^i(te+a) _^ gQ neglect all such oscillating terms (a rotating wave approximation). 
Then (jH^ reduces to a very simple matrix equation 



(36) 
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V V3{t) ] 
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The solution is easily obtained to be 



exp < 



i{V6-V2)hit 
24 























g-*01 
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cos(at) ze*'^isin(at) 
ie-'^'sm{at) cos(at) 



/J 















1 



MO) 

^i(O) 

MO) 
V ^3(0) J 



' MO) ^ 
MO) 
MO) 
\ MO) J 



^U{t) 



(37) 



^i(O) 
MO) 

\ MO) I 

where we have set a = ^^^ /j-i. That is, we obtained the unitary operator U{t). In 
particular, if we choose to and 0i satisfying 



cos(Q;to) = (sin(ato) = 1) and e*"^^ = i 



then 



-1 

1 



V 



(38) 



From this we want to construct the controlled NOT operator. However, it is not easy ^. 

At this stage we use a very skillful method due to Dirac [TU] . That is, we exchange 
two atoms in the cavity 



^U{to) is imprimitive in the sense of [S], so the main theorem in it says that our system is universal 
(namely, we can construct any element in C/(4)). However, how to construct a unitary element explicitly 



is not given in |H] 
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Exchange 



which introduces the exchange (swap) operator 



1 

1 



Multiplying U{to) by P gives 



PU{to) 



-1 



and from this we obtain 



(l2®(7i)PC/(io)(l2«)f7i) 



V 



-1 



(39) 



(40) 



This is just the controlled cr^ operator. Prom this it is easy to construct the controlled 
NOT operator, namely 



/ 

Cnot = (I2 ® W)C^^{l2 (^W)^ 
where W is the Walsh-Hadamard operator given by 

w = ^\ ^ ^ I = w-\ 

V2 1 _i 



1 

1 



(41) 
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See for example |n|. As to a construction of W by making use of Rabi oscillations see 
Appendix or [TT] . 

Therefore our system is universal jS], (Oj- 
A comment is in order. 

(a) If we choose to in U (t) satisfying 



cos{ato) 



{sm{ato) = 0) 



then we obtain the matrix 
/ 1 



V 



1 



CT3 (g) CT3. 



(b) In place of the ansatz (jHlj) we can set for example 

/ 



^ a++(t) ^ 



a+_(t) 





® |0) 






a_+(t) 

V u y \^ a__(t) ) 

Then we can trace the same line shown in this section and obtain a unitary operator 
under some resonance condition like ^E^. This is a good exercise, so we leave it to the 
readers. 



3.2 Controlled- Controlled NOT Gate 

In this subsection we treat the case of three atoms (the system of three qubits). To 
perform a quantum computation we need a (rapid) construction of the the controUed- 
controlled NOT gate. To construct it we need a construction of the controlled NOT gates 
of three types like 
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C NOT 



C-NOT 



C-NOT 



The controUed-controUed NOT gate (operator) 



CCnot ■■ C^®C^®C^ — >C^®C^®C^ 



is shown clS Qb picture 



CC 



NOT 



or in a matrix form 



V 



1 

1 



The (usual) construction by making use of controlled NOT or controlled U gates is shown 
as a picture ([0], [Sj) 
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where is a unitary matrix given by 



V 



I ( 1 + i 1-i 
1-i 1+i 



fo 1 
U 



(Ti. 



® |0). 



(42) 



However, we have not seen "reahstic" constructions in any references, so we give the 
exphcit construction. See the figure 2. 

To embed three-qubits in our quantized system we consider the simplest one as a wave 
function 

a+++{t) 
a++_(t) 
a+_+(t) 
a+__(t) 
a_++(t) 
a_+_(t) 

\ a— (t) J 

like in the case of two-qubits (jHT|) . 

To determine a dynamics that the coefficients a+++, a++_ 
tute into the equation 



, a satisfy we substi- 
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C-NOT 




Figure 2: The Controlled NOT gates of three types ((A), (B), (C) from the above) for 
the three atoms in the cavity 



22 



/ „itgBi/2 



V 



T^V{t)T 



J 



Let us rewrite the above equation. If we set 



\ 



-itgBj/2 



\ 



J 





/ \ 




















® |0) = 


























then we have 




^ gii3-Bi/2 









x/2 




a 


-+.{t)-a 






V2 






{t)+a+- + 


(t)-2a_++ 




Ve 




2a+_ 


_(i)-a_ + 


_(t)-a__+ 












f 


0.++- 


-(t)+a+_4 


.{t)+a-++ 










-(*)+«- + - 


-{*)+"— + 



v/3 

a___(t) 



|0) (43) 



■ d , , 



v 



/ 



v 



Q-itgBi/2 



Q-itgB^/2 



I 



(44) 



On the other hand, we have calculated e~**^^i/2 in ((T3I), e"'*3^3/2 (jjij) and the middle 
term Ttt>(t)T in 
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Therefore the (full) equation is reduced to the equations of {(fo, (pi, - ■ • , (p^} at the ground 
state |0) : 

/ n . . . \ / \ 



where 



dt 



ip4{t) 



Xi2 Xii Xi5 Xyj 

X21 X23 X2Q X2S 

X32 X34 xas X37 

2:41 X43 X410 a;48 

2^51 X53 X56 

Xq2 X64 2^65 ^67 

x-ji X73 0:76 X78 

a;82 a;84 ^87 



Mi) 
Mi) 

Mt) 
Mt) 



xu=pC{0)C{l), 
q — r 



g — r 

2^21 = Xi2, 
^23 — ^32, 
3^26 — ^62, 

q — r 



C(0)C(1), 



(/o(0)C(l) + 3Fo(0)5(l)), 



2^28 
2^34 



/-i(-l)C(0), 
C(0)C(1), 

mm, 



V2 

q — r 

—p + 2q + 2r 



X35 — X53, 

2p — q — r 



X37 = 



3^/2 



(/o(0)C(l) + 3Fo(0)5(l)), 



2^41 = 2^14, 
X43 — ^34, 
X4Q = Xq4, 



(45) 
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^48 = ^^~^~' /-l(-l)C(0), 

^5i = ^^(C(l)/2(2) + 5(l)Fi(2)), 

^53 = {C{1)M2) + 5(1)F,(2)) , 

2:56 = (/i(l)/2(2) + 4i/i(l)Fi(2) + 24/ii(l)/ii(2)) , 

^62 = ^^C(0)/i(l), 

-2^ + ^ + ^-C(0)A(i), 



(/o(0)/i(l) + 3Fo(0)i/i(l)), 



Xq4 


/o 

ov2 


^65 


— ^56) 


X67 


2{p + q + r) 
3 


X71 


= ^17, 


X73 


= ^37, 


X76 


= ^67; 


X78 


p + g + r 
V3 ^~ 


X82 


= 2^28, 


X84 


= ^48, 


X87 


= ^78 



/_i(-l)/o(0), 



and CiO), C(l), 8(1), /-i(-l), /o(0), /2(2), Fo(0), Fi(2), h{2), H,{1) are 

respectively given as 

p{t) = /i^e^^("i+'"^*+^i>, = /i2e'^("2+'"^*+^2>, r(t) = /i3e^^(^3+u;)t+03}^ 

C(0) = 1, 
C(l) = cos{tg), 
S{l)^sm{tg), 

/-i(-l) = 1, 

/o(0) = cos{tgV3), 
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„,^, 2cos(i^VlO) + 3 

h[^) - r ' 



i^o(O) 



/2(2) 




^i(l) 



F^{2) 



/3(2) 



/3(1) 




(iWlO- V73) 



These equations are complicated enough. 

Next let us derive the three controlled NOT gates (A), (B), (C) in Figure 2 from 
the matrix equation above. For that we use a resonance condition and rotating wave 
approximation associated to it. 

Derivation of (A) 

We focus on the components x^q and xq^ = x^q in the matrix. /i(l)/2(2) and hi{l)hi{2) 
contain the term Q-'^tay^o+Vf^ coming from the Euler formula, and therefore x^q contain 
the oscillating term ^ 



We note that this term is not contained in other components in the matrix. Here we set 
the resonance condition 



and apply the rotating wave approximation associated to this. Then the above compli- 
cated matrix equation becomes a very simple one 

'^Wo use p{t) in the first atom, however it is of course possible to use q{t) (the second atom) or r{t) 
(the third atom) 



gi{(ni +uj)t+4>i}^-itgy/ 10+^73 



gi{ (Hi +ui-gy/ 10+V73)t+<t>i } 




(46) 
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dt 
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(47) 



The solution is easily obtained to be 



exp < 
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ie''^isin(at) 




















ie~*'^isin(Q!t) 


cos (at) 
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\ / 



/J 



^ ^o(O) ^ 
^i(O) 

</^2(0) 
<^3(0) 
<^4(0) 
<^5(0) 
V^6(0) 
(^7(0) 



v 



¥^0(0) 
^i(O) 

<^2(0) 
<^3(0) 
(^4(0) 
<^5(0) 
<^6(0) 
<^7(0) 
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= UA{t) 



^ ¥^0(0) ^ 

^i(O) 

<^2(0) 

(^4(0) 

<^5(0) 
<^6(0) 
^7(0) 



where we have set a = ^q^^^' " 

Here if we choose Ia as cos(Q;iA) = — 1, then 
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By multiplying I2 ® cxi ® I2 from both sides we hav€ 



UA{tA) = (I2 0-1 (g) 12)C/a(^a)(12 <8) (7i ® I2) = 
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Moreover, by multiplying I2 <S) W <^ I2 from both sides we finally obtain the case (A) 
in Figure 2 



[l2^W^ l2)UA{tA){l2 IV ® I2) 



V 



1 

1 



Cnot^U- (50) 



Derivation of (B) 

We focus on the components 0:37 and 2:73 = 0:37 in the matrix. /o(0)C(l) and Fo(0)S'(l) 
contain the term Q-^^ai^+Vs) coming from the Euler formula, and therefore 0:37 contain the 
oscillating term 

gi{(ni+w)t+0i}g-itg(l+V3) _ ^i{(ni+uj-g(l+V3))t+(pi} 

This term is not contained in other components in the matrix. Here we set the resonance 
condition 

ni + (jj - g{l + \/3) = (51) 

and apply the rotating wave approximation associated to this. Then the above compli- 
cated matrix equation becomes a very simple one 



d 
I— 
dt 



I 



^2(1 - V3)/ii 



12 



00 00000 

00 00000 

e^^i 

00 00000 

00 00000 

00 00000 

e-^"^! 

00 00000 



\ 



(P2{t) 



(52) 
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The solution is easily obtained to be 



U)n(t) 
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<^6(0) 
V97(0) 



' MO) ^ 

<^2(0) 
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(^4(0) 
¥^5(0) 

<^6(0) 

997(0) 



(53) 



where we have set /3 



V2{-1+V3)hi 
12 
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Here if we choose Ib as cos{ptB) — —1, then 
/ \ 



UBitn) = 
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By multiplying I2 (8) I2 ® <Ji from both sides we have 



UsitB) = (I2 I2 «) cri)UBitB){l2 I2 ® (Ji) = I2 



V 



Moreover, by multiplying I2 ® I2 <8) Vt^ from both sides we finally obtain the case 
in Figure 2 



(I2 (8) I2 W)UB{tB){l2 ® I2 W^) = I2 



V 



1 

1 



— 12<S> Cnot- 



Derivation of (C) 

It is well-known that the case (C) is obtained as 



— 



® — ® 
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by making use of (A) and (B). However, 4 steps are needed in construction. We give 
a direct construction like (A) or (B) which is favorable from the point of view of quick 
construction. 

We focus on the components Xq^ and Xye = Xe7 in the matrix. /o(0)/i(l) and Fo(0)ifi(l) 
contain the term Q-^^giVs+VTo) coming from the Euler formula, and therefore Xgr contain 
the oscillating term 

gi{(ni+a;)t+0i}g-itg(V3+VlO) ^ ^i{(ni+u>-g{V3+VT0))t+,Pi} 

This term is not contained in other components in the matrix. Here we set the resonance 
condition 

+ cu - g{V3 + VlQ) = (55) 

and apply the rotating wave approximation associated to this. Then the above compli- 
cated matrix equation becomes a very simple one 
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(56) 
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The solution is easily obtained to be 



^ Mi) ^ 
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exp < 



it- 



-4 + V30)/ii 
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MO) 
MO) 
MO) 
MO) 
MO) 
MO) 



(57) 



where we have set 7 



(-4+V30)fti 
60 



33 



Here if we choose tc as cos{'-ytc 

I 



Uc{t 



c) 



) = —1, then 

1 

1 

1 

1 

1 

00000 -1 

00000 -10 

1 



By multiplying UA{tA) in ()49|1 from the left hand side we have 

'^lOOOO 

1 

1 

1 

1 

00000-10 

1 





Ucitc) = UA{tA)Uc{tc) 



V 



1 



Moreover, by multiplying I2 
in Figure 2 



[I2 ® I2 ® W)Uc{tc){'^2 ®12®W) 



W from both sides we finally obtain the case (C) 



'^lOOOOOOO^ 
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NOT- 



(58) 



34 



4 Further Problem 



In this section we consider the case of four atoms (the system of four qubits) and present 
a problem on constructing the controlled-controUed-controlled NOT gate defined by 



CCCnot 



as a picture. This gate is usually constructed as 




where is a unitary matrix given by 

/ 



V 



V 



3+i 1-i 
4 4 

1-i 3+i 
4 4 



1 

1 



(Ti. 



By the way, we have given the exact form of evolution operator for the four atoms 
Tavis-Cummings model |7^, therefore we can in principle track the same line shown in 
this paper and it may be possible to get the controlled-controlled-controlled NOT gate 
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or controlled-controlled-controlled unitary gates directly. However, such a calculation for 
the case becomes very complicated because we must treat 16 x 16 matrices at each step 
of calculations. We will attempt it in the near future. 

It is not easy to extend our method to the general case because to obtain the exact 
form of evolution operator for the general Tavis-Cummings model (see ©) is almost 
impossible. 

By the way, in [Hj we have given an idea on construction of the controlled NOT gate 
of type (A) by making use of the following skillful method (see Figure 3). 

(i) We move the third atom from the cavity. 

(ii) We insert a photon in the cavity as two atoms interact with it and subject laser fields 
to the atoms, and next exchange the two atoms, which gives the controlled NOT operator 
as shown in the preceding section. 

(iii) We return the third atom (outside the cavity) to the former position. 
It is easy to generalize our idea to the case of n atoms. 

One of important features of our model is that n atoms are trapped linearly in the 
cavity, which has both strong points and weak points. One of weak points is that a photon 
inserted in the cavity interacts with all atoms, so it is impossible to select atoms which 
interact with the photon. To improve this we presented the idea above.® 

To perform a full quantum computation we need to construct (many) controlled- 
controUed NOT gates or controlled-controUed unitary ones for three atoms among n- 
atoms, see jSj; §7. We believe that our idea is still available ^. 

In principle, we can construct general quantum networks. 

''We must estimate an influence of the "getting atoms (wliicli are not our target) in and out" on the 
whole states space, which is however difficult in our model. For that we must add in further terms 
necessary to calculate it 



36 




Figure 3: The process to construct the controlled NOT gate between the first atom and 
second one for the three atoms in a cavity 
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5 Discussion 



In this paper we constructed the controlled NOT operator in the system of two qubits 
and controUed-controlled NOT operator in the system of three qubits in the quantum 
computation based on Cavity QED, which show that our system is universal. Therefore we 
can in principle perform a quantum computation. Our method is completely mathematical 
and we use several Rabi oscillations in a consistent manner. 

We hope that some experimentalists will check whether our method works 
good or not. 

See [IF and their references for some experiments on Cavity QED (which may be related 
to our method). 

We conclude this paper by making a comment (which is important at least to us). The 
Tavis-Cummings model is based on (only) two energy levels of atoms. However, an atom 
has in general infinitely many energy levels, so it might be natural to use this possibility. 
We are also studying a quantum computation based on multi-level systems of atoms (a 
qudit theory) ^21- Therefore we would like to extend the Tavis-Cummings model based 
on two-levels to a model based on multi-levels. This is a very challenging task. 

Acknowledgment. We wish to thank Tatsuo Suzuki and Shin'ichi Nojiri for their helpful 
comments and suggestions. 

Appendix 

One Qubit Operators by Classical Fields 

Let us make a brief review of theory without the radiation field, whose states space is 
only tensor product of two level systems of each atom. See the figure 5. 
The Hamiltonian in this case is 

^ = E ^ + (aj+)e*(^^*+<^^) + aj-)e-^(^^*+^^)) } (59) 
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Figure 4: The n atoms in the cavity without a photon (in Figure 1) 

from (0), where we have omitted the unit operator 1 for simphcity. Here let us remind 
the notation 



Mj = I2 



Then 



A 

2 



I2 ® M ® I2 ® ■ ■ • ® I2 for l<]<n. 



E 



I e 2 



e 2 



^ f h, \ I e-^ 



e 2 



/J 



(f/i ® ■ ■ • ® f/„) 5] 

.=1 \ /i, 



2 J 



where 



A 
2 



(60) 



e 2 



V 



e~* 2 



The wave function defined by = H\^!) with can be written as a tensor 

product 

1^) = l^l)®---® l^n), (61) 

SO if we define 
then it is easy to see 
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A-n 



Th solution is easy to obtain 



=(g)exp<( -zt 



2 

hi 



V 

Therefore, the solution that we are looking for is 



1^(0)). 



' e 2 



.=1 V 

Last we note that 



e~* 2 



exp < —it 



exp < — it 



A-f7, 



1^(0)). (62) 



2 



y X21 X22 



where 



Xn = cos \ 't\ — + h'^ 
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i— 
2 



sin ( tx/f + 



'i + h^ 



X12 = X21 = -ih- 



sin U '4 + 



'i + h^ 



X22 = COS \t\ — + \ + t- 



sin 



02 



We can always construct unitary operators in U (2) at each atom by using Rabi osillations, 
see for example jTTj . 

As an example let us construct the Walsh-Hadamard operator W which has been used 
in the text. We set 



Vit) 



I e 2 



V 



_ ■ nt+4> 

e 2 
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j exp < 


—it 

. V 



/ A- 



h 



\ 



h - 



A-n 



and set the resonance condition A = Q, then 
V{t) 



I e 2 



V 



e 2 



e 2 



/ 



-j(At+0) 



cos{ht) —isin{ht) 
-isin{ht) cos{ht) 

cos{ht) —ism{ht) 
—ism{ht) cos{ht) 
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Now we again set 



Vit) 



V 



-i(Ai+(/)) 



^ cos{ht) —ism{ht) 
^ ~ism{ht) cos{ht) 



(63) 



for simplicity because there is no interest in the overall phase factor ^. On the other hand, 
when we don't subject a laser field to the atom {h = 0) we have 



Vo{t) 



lAt 



(64) 



By using V{t) and Vo(t) we construct W in the following ([H]). First we set 

V{tf,U) = Vitf-U), Voitf,ti) = Voitf-U) for tf > U. 
A sequence of operators to construct the Walsh-Hadamard gate is given as follows 



Vn 



V 



Vo 



ti 



t2 



time 



For Vo(ti,0), Vo(t3,t2) with ti = 3n/2A and 1^-1^ = 37r/2A 



^(ti,0) = 
and V{t2.,ti) with t2 — ti = n/Ah 



1 
i 



-i{A{t2~ti)+(f 



1 
i 



V 73 72 



we have 



Voit3,t2)Vit2,h)Voiti,0) 



V 



1 


1 


V2 




1 


1 


V2 


V2 



^In the Hamiltonian a constant term which makes a overaU phase factor has been removed 
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By choosing the phase as e ^■t^^*^ ti)+<i>} — i finally obtain 

Vo{h,h)V{t2,h)Vo{h,0) = -^\ ^\=W. (65) 
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